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We present a quantitative finite temperature analysis of a recent experiment with Bose-Fermi mixtures in op- 
tical lattices, in which the dependence of the coherence of bosons on the inter-species interaction was analyzed. 
Our theory reproduces the characteristics of this dependence and suggests that intrinsic temperature effects play 
an important role in these systems. Namely, under the assumption that the ramping up of the optical lattice is an 
isentropic process, adiabatic temperature changes of the mixture occur that depend on the interaction between 
bosons and fermions. Matching the entropy of two regimes — no lattice on the one hand and deep lattices on the 
other — allows us to compute the temperature in the lattice and the visibility of the quasi-momentum distribution 
of the bosonic atoms, which we compare to the experiment. We briefly comment on the remaining discrepancy 
between theory and experiment, speculating that it may in part be attributed to the dependence of the Bose-Fermi 
scattering length on the confinement of the atoms. 



I. INTRODUCTION 

Ultracold atoms in optical lattices are, due to the available 
impressive control over system parameters, ideal candidates 
for "quantum simulators" that mimic condensed matter sys- 
tems d [2). We have already seen them display the transi- 
tion from a superfluid to a Mott insulator [3], Fermi surfaces 
have been observed [4], and, recently, the finite temperature 
phase diagram for bosonic superfluids in an optical lattice has 
been obtained experimentally \5 \. Multi-component mixtures, 
among them mixtures of bosonic and fermionic atoms, offer 
a variety of additional quantum phases of matter. Charge- 
density waves 0, supersolids 00, exotic superfluid l9l[T0l 
and Mott-insulator phases fTTlfT3ll have been predicted and 
we will certainly see experimental signatures of these in the 
near future. Of course, temperature plays a prime role for such 
quantum simulators and its influence needs to be understood 
or, better yet, be under control. But, even just determining 
the temperature in an optical lattice is an extremely difficult 
task Q. Thermometry methods for such systems without lat- 
tices are however well established. Hence, under the assump- 
tion that the lattice is ramped up adiabatically (which is usu- 
ally a good approximation and was also recently confirmed 
for bosons in optical lattices |5|), i.e., without changing the 
entropy, one may infer about the temperature in the lattice us- 
ing entropy-matching methods. Not only does this hold the 
opportunity for thermometry in the lattice but also offers the 
possibility to further cool the atoms |[T4l . 

In this work, we study interaction-dependent temperature 
effects occurring in Bose-Fermi mixtures under this assump- 
tion of raising the optical lattice being isentropic. We compare 
our results to the visibility of the quasi-momentum distribu- 
tion, which was obtained in the experiment in Ref. [15 ] for a 
87 Rb- 40 K mixture. By matching the entropy of two very dif- 
ferent regimes (with and without lattice), we are able to take 
all experimental parameters (such as the anisotropic trapping 
potential, the number of particles and the lattice parameters) 
into account, leaving no free parameters in our theory. Our re- 
sults show that these effects have a significant influence on the 
coherence of the bosonic atoms and depend quite strongly on 
the interaction between the two species, in agreement with the 
experiment. Hence, we are faced with a situation in which in- 




FIG. 1 : Entropy as a function of temperature and Bose-Fermi scat- 
tering length clbf (oc qbf oc V, the Bose-Fermi interaction) in units 
of the Bohr radius ao with (translucent surface) and without (opaque 
surface) a twelve recoil energies deep optical lattice. Bold lines at 
clbf — —400 show the same for a purely bosonic system, for which 
the entropies are of course independent of clef (upper line: includ- 
ing the lattice, lower line: without lattice). The mixture consists of 
4 x 10 5 87 Rb and 3 x 10 5 40 K atoms and trap parameters are as in 
the experiment in Ref. [15], see footnote 1 16 1 for details. 



trinsic adiabatic temperature effects play a dominant role and, 
as we argue below, have already been observed experimen- 
tally. 

More specifically, we calculate the entropy as a function 
of the temperature in the absence of the lattice invoking 
the Hartree-Fock-Bogoliubov-Popov mean-field approxima- 
tion for the bosons, including the fermions in a self-consistent 
mean field interaction. To describe the mixture in the lattice, 
we use the single-band Bose-Fermi-Hubbard Hamiltonian and 
calculate the entropy as a function of the temperature for a 
deep lattice perturbatively. For both regimes, we assume that 
the mixture is in thermal equilibrium such that we can assign 
one temperature to it. This results in temperature-entropy di- 
agrams as in Fig. [T] and allows us to obtain the temperature 
in the lattice Tf as a function of the initial temperature Ti by 
matching the corresponding entropies, see Fig. [2] The result- 
ing adiabatic heating or cooling of these isentropic processes 
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was already analyzed for purely fermionic ifTTl and bosonic 
lfT4l [T8l systems, a Fermi gas of atoms that can pair into 
molecules via a Feshbach resonance lfT9l , and Bose-Fermi 
mixtures | 20 ] as in this work. For fixed inter-particle interac- 
tion, loss of bosonic coherence due to presence of fermions 
was observed in the experiments in Refs. I2T1 [22l and at- 
tributed to intrinsic temperature effects in Refs. I2Ql [2H |23) 
(see, however, also Refs. (2H [23, in which different expla- 
nations were put forward). The adiabatic assumption was re- 
cently confirmed experimentally for bosons [5 ]. These studies 
show that Tf can depend strongly on the system parameters, 
prime among them the inhomogeneity introduced by the trap- 
ping potential, and hence realistic descriptions should take 
all of them into account. To the best of our knowledge, the 
present analysis is the first to be directly comparable to the ex- 
periment over the full range of inter-species interactions and 
to fully take the experimental situation into account. 

As our scheme to obtain Tf is similar to the one used in 
Ref. ll20l (note however, that here we do not take geomet- 
rical averages of the trapping frequencies but take the full 
anisotropic situation into account and do not invoke the lo- 
cal density approximation in the lattice), we merely outline it 
in the two subsequent Sections [IT] and III Being equipped with 
Tf, we calculate the visibility of the quasi-momentum distri- 
bution within thermal perturbation theory and compare it to 
the experiment [15 ] in Section [V] resulting in Fig. [3] 

II. NO LATTICE 

We start from the microscopic model of a mixture of 
bosonic atoms of mass tub and fermionic atoms of mass 
tuf subject to respective trapping potentials Ve(r), Vp(r), 
and interacting via contact interactions parametrized by the s- 
wave scattering lengths clbb (Bose-Bose interaction) and cibf 
(Bose-Fermi interaction), 

H=J2 f*r%(r)[=!£p-+V s (r)-» 8 ]* s (r) 

- ^Jdr^ B (r)^ B (r)^ B (r)^ B (r) W 
• g B F \ dr *^(r)* B (r)*^(r)* F (r). 



S=B,F^ 



Invoking the Hartree-Fock-Bogoliubov-Popov mean- field ap- 
proximation (26l |23 for the bosons and including the 
fermions in a self-consistent mean-field approximation (23 
[22), yields the following set of coupled equations, which 
we solve self-consistently with an iterative numerical scheme 
(see, e.g., Ref. |20] for more details), (i) The finite tempera- 
ture Gross-Pitaevskii equation in the Thomas-Fermi approxi- 
mation for the condensed atoms, 



n (r) 



max 



{0, 



HB-VB{r)-g BF m(r) 
9bb 



2n T (r)}, (2) 



where m(r) is the fermionic density, tit(t) the den- 
sity of thermal bosons, qbb = 47rft 2 a^/m^, qbf = 
2irh 2 (iBF(mB + trf) / (mBrriF), and the chemical potential 



Hb controls the total number of bosons Nb = Jdrn(r), 
n(r) = no(r) + riT(r). (ii) The thermal contribution 



, (3) 



where /3 = l/^sT) denotes the inverse temperature and the 
Bogoliubov amplitudes and quasi-particle spectrum are given 
by 

/ \ t£^+V B (r)-/j, B +2g BB n(r)+g FB rn(r) -, 
u ±fa r ) = — 2^r) ± 2> 

*(P> r ) = E + V B (r) - m + 2g BB n(r) ( 4 > 
+ gFBm(r)] 2 -g 2 BB nl(r). 

(iii) The fermionic density in local density approximation 



m(r) 



dp 



(2tt) 3 e ^S( P , r ) +1 i 



where 



(5) 



(6) 



and hf controls the total number of fermions Nf = 
J drra(r). 

After convergence, we are in a position to compute the 
entropy of the mixture in thermal equilibrium, S(T)/kB — 
j dpdr[s B (p : r) + r)]/(27r) 3 , with individual contri- 

butions given by 



sf(p,v) - e ^ (p ,, )+1 



e Pe(p,r)_l 



log (1 
log (1 



(7) 



For the parameters fT6l of the experiment in Ref. fT5lL we 
show the resulting entropy as a function of the temperature in 
Fig. [T] We can see that it depends only weakly on the inter- 
action gBF - in stark contrast to the situation including the 
lattice, which we treat in the subsequent Section III Further- 
more, over the whole range of interactions, it is higher than 
the entropy for a purely bosonic system (note that this is not 
the same as the non-interacting situation, simply due to the 
fermionic contribution sf to the entropy of the mixture), for 
which the entropy may be obtained from the approximative 
expression for the energy l26l 



N B k B T° 



at* 



b[N B /6 ^(l-t 3 )] z/0 (5 + l6t 3 ), (8) 



2/5/ 



which holds over a wide temperature and parameter range 
E81I26II23. Here > a = 3C(4)/C(3), b = C(3) 1/3 15 2 / 5 /14, 
aho is the harmonic oscillator length of the trapping poten- 
tial, and T® the critical temperature for bose condensation in 
the absence of interactions. This nourishes the hope of being 
able to also arrive at a closed expression for the entropy of the 
mixture, which is however beyond the scope of this work. 
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III. DEEP OPTICAL LATTICE 

For sufficiently deep optical lattices, the system may be de- 
scribed l30l by the single-band Bose-Fermi-Hubbard Hamil- 
tonian ED, H = J + H , 

J=-JB^b\b j -J F 'Z f f}f j , 

(iJ) (h3) (9) 

Ho = ^ [Uhi(hi — l)/2 + Vhifhi — Hifii — Viihi] , 

i 

where hi = h\hi, rhi = fjfi, the hi (fi) are bosonic 
(fermionic) annihilation operators, and the site-dependent 
chemical potentials, Hi = n — tf , Vi = v — tf , account for 

the trapping potentials tf^ F and control the number of parti- 
cles. This model is obtained by including the lattice into the 
microscopic model in Eq. ([I]), expanding the field operators 
in the Wannier basis, and neglecting all bands above the low- 
est band [30] and contributions beyond nearest neighbours. 

B IF 

The amplitudes Jbif^^i ^ oc 9bf oc clef are then 
obtained from a single-particle band- structure calculation for 
appropriate lattice parameters fTGl and the overlap of the re- 
sulting Wannier functions. For a discussion of the validity of 
the single band approximation and contact interaction, we re- 
fer the reader to Refs. l24l [32ti35l 

Similar to Ref. l20lL we obtain the entropy as a func- 
tion of temperature employing the thermodynamic interac- 
tion picture, treating J as a perturbation up to first order (the 
difference to Ref. [20] being that we do not employ a lo- 
cal density approximation and do not take geometrical av- 
erages of the trapping frequencies), which is applicable for 
1 > Jb/fP (- nK/T in the situation at hand) [36]. This 
amounts to approximating e -/3 ^ w e _/3 ^°(l — 7), 7 = 
Jq dx e x ^° Je~ x ^ Q . For given temperature this yields the par- 
tition function Z = tr[e~^^] « tr[e _/ ^° (1 - 7)] and the total 
number of bosons and fermions as a function of the chemical 
potentials. Numerically solving (J^r^) = Nb, rhi) = 
N F for /i, v then yields the entropy S(T) = /3(H) + log(Z) 
for given temperature and particle numbers N B / F up to first 
order in J. Note that we consider the full three-dimensional 
anisotropic experimental situation, i.e., we need to take a large 
number of lattice sites into account, which results in quite in- 
volved numerics. 

Fig. [T] summarizes the result of this procedure for the ex- 
perimental parameters of Ref. |[T5l . The most prominent fea- 
ture of S(T) is its strong dependence on the Bose-Fermi in- 
teraction V oc aBF, while we found only a weak dependence 
in Section [TT| For fixed temperature, starting from a plateau 
for strong attraction, the entropy increases until it reaches a 
maximum around aBF = 0, from which it decreases with 
increasing a^F to a plateau for strong repulsion. It is this be- 
haviour that will crucially influence the temperature Tf in the 
lattice and hence also the coherence properties of the bosonic 
atoms, which displays the same strong dependence on the 
Bose-Fermi interaction, see Section [V] The plateaus for large 
\a BF \ are easily explained: For large repulsion, phase sepa- 




FIG. 2: Temperature of the mixture in an optical lattice (obtained 
by entropy matching from the data in Fig. [TJ as a function of the 
Bose-Fermi scattering lenght aBF (oc qbf oc V, the Bose-Fermi 
interaction) and the initial temperature without lattice . The bold 
line at aBF = 400 shows the same for a purely bosonic system, for 
which the final temperature is of course independent of as f • 



ration occurs and once this phase is entered, increasing a B F 
further does not have any effect. For large attraction on the 
other hand, bosons and fermions are forced to occupy the 
same lattice sites, and again further increasing \clbf\ does 
not have any effect. Comparing the entropy of the mixture 
to the purely bosonic situation (note that this is not the same 
as the non-interacting case as, of course, also for a^F = 0, 
the fermions contribute to the total entropy), we see that the 
former is always higher than the latter for the considered pa- 
rameter regime. As we will see in the subsequent Section, 
while the adiabatic ramping up of the lattice leads to adiabatic 
cooling, this causes the mixture to be less cooled than bosons 
would be without fermions. 



IV. THE TEMPERATURE IN THE LATTICE 

Having obtained the entropies with, Sf, and without, Si, 
lattice, we are now in the position to obtain the temperature 
in the lattice, Tf, for given initial temperature, Ti, by match- 
ing the respective entropies Sf(Tf) = Si(Ti). If the opti- 
cal lattice is indeed raised adiabatically and the mixture is in 
thermal equilibrium, this enables us to compute Tf as a func- 
tion of Ti, which can be measured as, without lattice, ther- 
mometry methods are well established. Fig. [2] shows the re- 
sult obtained by matching the entropies in Fig. [T] As Sf in 
Fig. [T] already suggests, we find a strong dependence of the 
temperature in the lattice on the Bose-Fermi interaction, most 
pronounced for high initial temperatures. The qualitative be- 
haviour of Tf is similar to Sf: For fixed Ti starting from a 
plateau at large attraction, the temperature decreases with in- 
creasing aBF, reaches a minimum around aBF = 0, from 
which it increases with increasing aBF to a plateau at high re- 
pulsion. We also depict Tf for a purely bosonic system, which 
shows that while over most of the parameter regime, the mix- 
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FIG. 3: Visibility, V, of the time-of-flight distribution at a lattice 
depth of twelve recoil energies as a function of the Bose-Fermi scat- 
tering length clef (oc qbf oc V, the Bose-Fermi interaction) and 
initial temperature without the lattice T\ (left surface, Tf obtained 
as in Fig. [2]> and temperature in the lattice Tf (right surface). Inset 
shows the same at the initial temperature indicated by the bold line 
on the left surface for a purely bosonic system (including an uncer- 
tainty of 10% in the initial temperature) and the mixture (including 
an uncertainty of 15% in the initial temperature). The initial temper- 
ature was chosen such that the purely bosonic situation matches the 
experiment 1 15|. Circles are the data points from said experiment. 
The mixture consists of 4 x 10 5 87 Rb and 3 x 10 5 40 K atoms and 
all parameters (trap geometry, lattice depth) are as in this experiment 
and there are no free parameters in the theory. 

ture is cooled, the cooling is less than it would be without 
fermions. Being equipped with Tf, we can now move on to 
study the dependence of the bosonic coherence on the Bose- 
Fermi interaction. 



V. VISIBILITY OF THE BOSONIC QUASI-MOMENTUM 
DISTRIBUTION 

After switching off all potentials and letting the atom cloud 
evolve freely for a time t, the density of bosons is well ap- 
proximated bv l37ll38l 

n(p) oc ^(SJS^e^-iJei^dil 2 -!*! 2 ) 

ij (10) 

where r = Htj '(2ra#a 2 ) and w is the Fourier transform of the 
Wannier function centered at zero. This density is measured 
by taking an absorption image of the cloud, resulting in the 
column density n(p x ,p y ) = J dp z n(p). For shallow lattices 
and low temperatures, i.e., in the superfluid regime, this den- 
sity displays a pronounced interference pattern, which van- 
ishes deep in the Mott regime (ultra-deep lattices) and for high 
temperatures. Hence, the visibility of this interference pattern, 
V = (n max - n min )/ (n max + n min ) 1 39 ], is an indicator for 
the coherence of the 87 Rb atoms. We calculate V for a given 
temperature and number of particles by computing the two- 
point correlations (b\bj) in thermal perturbation theory up to 
first order in J, (b\bj) « tr[e-^°(l -7)&t^]/tr[e-^°(l - 



7)] . From this we obtain the column density n(p x , p y ) and the 
visibility V |39) up to first order in J. The result of this com- 
putation is shown in Fig. [3] for parameters as in Ref. lfT5lfT6l 
and for two different temperatures: We depict V as a function 
of the temperature in the lattice Tf (right surface) and as a 
function of the initial temperature Ti (left surface). This cor- 
responds to two different scenarios: (a) A scenario in which 
adiabatic cooling mechanism are omitted and hence the tem- 
perature in the lattice Tf does not depend on the Bose-Fermi 
interaction and (b) the more realistic scenario in which the 
final temperature does depend on through the entropy 
matching described above. As Fig.[3]shows, the two scenarios 
result in opposed behaviours of V as a function of the interac- 
tion strength clef- For all temperatures and large \clbf\, the 
visibility V is higher on the attractive side than on the repul- 
sive side of the interaction for scenario (a), while we see the 
exact opposite for scenario (b). 



A. Comparison to the experiment 

Before we compare our results to the experiment, we re- 
call the assumptions and approximations of our theory. We 
have computed the temperature in the lattice by matching the 
entropies of the mixture with and without lattice (under the 
isentropic assumption, i.e., under the assumption that the lat- 
tice is ramped up adiabatically, which was recently confirmed 
for 87 Rb atoms |5], this indeed yields the temperature in the 
lattice as a function of the initial temperature). To obtain the 
latter, we invoked the Hartree-Fock-Bogoliubov-Popov mean- 
field approximation, and computed the former within thermal 
perturbation theory up to first order in the tunnelling param- 
eter of the single-band Bose-Fermi-Hubbard model (note that 
while the influence of higher bands in the lattice can not be 
completely ruled out and can have an effect on the bosonic 
coherence l24l [32), we are in a regime in which their occu- 
pation is expected to be small (30)). For all calculations, we 
assumed the mixture to be in thermal equilibrium such that we 
can assign one temperature to the mixture. This assumption 
is well justified if the mixture has sufficient time to equili- 
brate and gets worse the smaller the Bose-Fermi interaction. 
Finally, we computed the visibility of the time-of-flight inter- 
ference pattern, within thermal perturbation theory up to first 
order in the tunnelling parameter. 

The inset of Fig. [3] compares our results directly to the ex- 
periment in Ref. IU5F . We assume that the initial temperature 
for all measurements was approximately the same as in the 
purely bosonic situation. We can see that the results display 
the same qualitative behaviour: Starting at strong attraction, 
the visibility decreases to a minimum, from which it increases 
with increasing up to a maximum at around clef = 0, 
and finally decreases to a plateaux for strong repulsion. This 
is in stark contrast to what one finds without taking intrinsic 
temperature effects into account (see right surface in Fig. [3}: 
The visibility would, at the relevant temperatures, simply de- 
crease monotonically with increasing interaction. 

Fig. [3] also shows a shift of the theoretical results relative 
to the experimental data. All our results are parameterized by 
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the Bose-Fermi scattering length clef* which is tuned in the 
experiment by addressing the magnetic Feshbach resonance at 
around B w 546.9 G (40} [4T]| . Close to resonance, magnetic 
field B and scattering length are related by 

a BF = af F (1 - (11) 

where the resonance is at Bo, a^p * s me background scatter- 
ing length, and AB the width of the resonance. A faithful 
experimental determination of all the parameters in this re- 
lation is an extremely difficult endeavour. They depend on 
external parameters such as the trapping potential l34l - and, 
even more so, the tight "trapping" within a lattice site l33ll35l . 
Especially the latter complicates a direct comparison to the 
experiment: Neither ab initio calculations for realistic inter- 
atomic potentials nor experimental measurements of the de- 
pendence of aBF on the magnetic field are available for the 
situation at hand. In fact, the experimental visibility shown in 
Fig.[3]is really a function of the magnetic field, with clef de- 
termined from Eq. {IT) with B = 546.9 G, AB = -2.9 G, 
and a bp = — 185ao (42). Of course, due to the mentioned 
difficulties, this scattering length is a priori not the same as 
the one used to model the inter-atomic contact potential. Ex- 
ploring wether this may explain the remaining discrepancy be- 
tween theory and experiment in Fig. [3] is an exciting theoret- 
ical and experimental challenge that might result in a better 
understanding of the dependence of aBF on external poten- 
tials. We hope that, due to the pronounced features of the 
visibility — most prominently the location of the maximum — 
the present study can contribute to the work along this direc- 



tion. 



VI. CONCLUSION 

We have studied intrinsic temperature effects in Bose-Fermi 
mixtures taking the full three-dimensional anisotropic exper- 
imental situation into account. Under the adiabatic assump- 
tion, we have determined the temperature in the lattice as a 
function of the temperature before the lattice is ramped up 
and found a strong dependence on the inter- species interac- 
tion. This dependence affects the coherence of the bosons 
and is displayed in the visibility of the time-of-flight interfer- 
ence pattern, which we have compared to the experiment in 
Ref. 1 15 ], finding qualitative agreement. Not including these 
temperature effects results in a very different dependence on 
the interaction and lead us to conclude that they need to be 
incorporated into any realistic description of Bose-Fermi mix- 
tures in optical lattices. The remaining discrepancy between 
theory and experiment could be attributed to the dependence 
of the inter- species scattering length on external (trapping) po- 
tentials and we hope that it motivates further investigations of 
said dependence. 
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